Explosive transitions to synchronization in networked phase oscillators 
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We introduce a condition for an ensemble of networked phase oscillators to feature an abrupt, 
first-order phase transition from an unsynchronized to a synchronized state. This condition is met 
in a very wide spectrum of situations, and for various oscillators' initial frequency distributions. We 
show that the occurrence of such transitions is always accompanied by the spontaneous emergence 
of frequency-degree correlations in random network architectures. We also discuss ways to relax the 
condition, and to further extend the possibility for the first-order transition to occur, and illustrate 
how to engineer magnetic-like states of synchronization. Our findings thus indicate how to search 
for abrupt transitions in real- world applications. 
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Critical phenomena in complex networks, and the 
emergence of abrupt dynamical transitions in the macro- 
scopic state of a system are currently a subject of the 
utmost interest [1-4 . As far as the synchronization of 
an ensemble of networked phase oscillators is concerned 
[5Hll , the occurrence of a first-order phase transition was 
described, so far, in two rather special and apparently 
opposite situations: an all-to-all coupling architec- 
ture, with oscillators' frequencies evenly spaced, that is, 
a measure-zero realization of a uniform distribution [8], 
and a) scale- free connection topology, where a positive 
correlation between the heterogeneity of the connections 
and the frequencies of the oscillators is introduced ad hoc 

mm- 

In this Letter, we show how a sharp, discontinuous 
phase transition can occur, instead, as a generic feature 
of the synchronization of networked phase oscillators. 
Precisely, we initially give a condition for the transition 
from unsynchronized to synchronized states to be first- 
order, and demonstrate how such a condition is easy to 
attain for any oscillators' initial frequency distribution. 
We then show how such transitions are always accompa- 
nied by the spontaneous emergence of frequency-degree 
correlation features. Third, we show that the considered 
condition can be even softened in several cases. Finally, 
we illustrate, as a possible application, the option of ex- 
pressing magnetic-like states of synchronization with the 
use of such transitions. 

We consider a network of Kuramoto [11] oscillators: 

d(j)- ^ 

= uji^ d^^aij^\n{(t)j -(t)i), (1) 

where (j)i is the phase of the i^^ oscillator {i = 1, A^), uJi 
is its associated natural frequency drawn from a generic 
frequency distribution p(cj), d is the coupling strength, 
and {aij} are the elements of the adjacency matrix that 
uniquely defines the graph. The classical order parameter 
for system (jlj) is r{t) = | Xl^li e*"^^*^^^ |, and the level 



of phase synchronization can be monitored by looking at 
the value of S = {r{t))T^ where (. . . )t denotes a time 
average with T ^ 1. Furthermore, for each oscillator z, 
we denote by JV{i) the set of oscillators linked to it. 

As the coupling strength d increases, system ([T]) un- 
dergoes a phase transition from the unsynchronized {S ^ 
l/^/N) to a synchronous {S ^ 1) state, where all oscil- 
lators ultimately acquire the same frequency. For this 
phase transition to display a first-order feature, we have 
to avoid that any oscillator behaves as the core of a clus- 
tering process, where its neighbors begin to aggregate to 
the synchronous state smoothly [12 . This could be real- 
ized by considering a certain frequency difference with its 
neighbors. We have observed that a condition in order 
to achieve a first order phase transition is: 

(A) for each oscillator i, all nodes j belonging to A/'(i) 
satisfy \uji — ujj\ > 7c 

Condition (A) is tantamount to imposing a minimal value 
for the frequency difference between linked nodes. We 
now fix = 500, and illustrate the synchronization route 
for several frequency distributions p{oj) when condition 
(A) is met. To this purpose, we consider network topolo- 
gies resulting from the following procedure: i) we assign 
natural frequencies {a;^}, drawn from a distribution p(co'), 
to the N oscillators; ii) we randomly pick a pair 
of oscillators, and form a link between them only if the 
value of \uji — ujj\ exceeds a given gap 7; Hi) we repeat 
point ii) until the desired number of links L in the graph 
is formed. After a final check on the connectedness of 
the resulting network, the procedure yields Erdos-Renyi- 
like [13] topologies with an average degree (k) = 
We then use the obtained adjacency matrix to simulate 
system ([T]), and monitor the state of the network as a 
function of the coupling strength (i, by gradually increas- 
ing the value of d in steps 6d = 10""^, from d = 0. At 
each step, a long transient is discarded before the data 
are acquired for further processing. Moreover, insofar 
as we are looking for a first-order phase transition (and 
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FIG. 1: (Color online) (Top panels) Phase synchronization 
level S (see text for definition) vs. the coupling strength d, 
(left panel) for different values of the gap 7 and (k) = 40, 
and (right panel) different values of the average degree (k) 
with 7 = 0.4. The frequencies of the oscillators are here 
drawn from a uniform distribution within the range [0,1]. 
See text for the construction procedure of the networks. In 
both panels, the continuous (dashed) lines refer to the forward 
(backward) simulations. (Bottom panels) S in the parameter 
space (c?, 7). The values are color coded, according to the 
color bars. (Left panel) (k) — 20 and (right panel) {k) — 60. 
The horizontal dashed lines mark the separation between the 
region of the parameter space where a second-order transition 
occurs (below the line) and that in which the transition is 
instead of the first order type (above the line). The striped 
area delimits the hysteresis region. 



thus for an expected associated synchronization hystere- 
sis), simulations are also performed in the reverse way, 
i.e. starting from a given (imax (where the ensemble is 
phase synchronized), and gradually decreasing the cou- 
pling by 5d at each step. In what follows, the two sets 
of numerical trials are termed as forward and backward^ 
respectively. 

In Fig. [l]we report the results obtained by setting p{(jj) 
as a uniform frequency distribution in the interval [0, 1]. 
The top panels depict 5* as a function of d. Namely, 
the top left (right) panel illustrates the case of a fixed 
mean degree {k) = 40 (of a fixed frequency gap 7 = 
0.4), and reports the results of the forward and backward 
simulations at different values of 7 {{k)). The bottom 
panels show S in the parameter space d — 7 for (k) = 20 
and (k) = 60. 

Several pieces of information can be extracted from the 
Figure. First, the rise of a first order phase transition is 
a generic feature in the parameter space: the horizontal 
dashed lines in the bottom panels mark the value of 7c, 
separating the two regions where a second-order transi- 
tion (below the line) and a first-order transition (above 
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FIG. 2: (Color online) (Top panels) Degree ki that each node 
achieves after the network construction is completed vs. its 
natural frequency uJi. (k) = 100 and frequency gaps (a) 7 = 
0.0, and (b) 7 = 0.4. The red solid line in panel (b) represents 
the theoretical prediction /(cj) (see text). (Bottom panels) 
Average of the local natural frequency (uj) for j G A/'(i) vs. 
the natural frequency Ui of the i^^ node of the network for 
(c) 7 = 0.0 and (d) 7 = 0.4. 



the line) occurs. Second, while increasing (k) of the net- 
work facilitates the occurrence of the explosive transition, 
as both the values of 7 and d for which the first-order 
phase transition takes place shrink, it also reduces the 
region of hysteresis. This is consistent with the results of 
Ref. [8] for an all-to-all connected case, where in the limit 
^ 00 a first-order phase transition has been predicted 
in the absence of hysteresis. 

Another relevant result is the spontaneous emergence 
of degree-frequency correlation features associated to the 
passage from a second- to a first-order phase transition. 
While such a correlation was imposed ad hoc in Refs. 
[9l Ho], here condition (A) creates for each oscillator i 
a frequency range around a;^, where links are forbidden. 
The final degree ki is proportional to the total probability 
for that oscillator to receive connections from other oscil- 
lators in the network, and therefore to 1 — J^'^^ p{uj^)duj\ 
This is shown in Fig. |2|a)-(b), where the degree ki that 
each node achieved after the network construction is com- 
pleted is reported as a function of its natural frequency 
uji^ for (k) = 100. Panel (a) refers to the case 7 = 
in which no degree-frequency correlation is present in 
the resulting network. In panel (b), instead, we re- 
port the case 7 = 0.4 (a value for which a first-order 
phase transition occurs) and the (conveniently normal- 
ized) function f{uj) = 1 — j'^^'^ p{(jj')duj' ^ with p{(jj) = 1 
for uj e [0,1], and p{uj) = elsewhere, which gives ev- 
idence of the emergence of a very pronounced V-shape 
relationship between the frequency and the degree of the 
network's nodes. Inspecting the average frequency of 
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FIG. 3: (Color online) S vs. d resulting from the forward 
(continuous lines) and backward (dashed lines) simulations of 
system ([T]), displaying first-order phase transitions for differ- 
ent initial frequency distributions ^(cj), defined on the interval 
[0,1] (left panels), and distribution of the final node degree 
ki vs. the corresponding oscillator's natural frequency (right 
panels). From top to bottom: (a)-(b) Rayleigh distribution 
for 7 = 0.3. In (b), the red solid line represents the theo- 
retical prediction /(cj); (c)-(d) uniform distribution using the 
local mean field condition for two values of 7. In panel (d) 
7 = 0.4; (e)-(f) Gaussian distribution with Z = 0.7; (g)-(h) 
bimodal Gaussian with Z = 0.8. The insets in panels (a), (e), 
and (g) report the corresponding distribution p{uj). See text 
for the details on the specific construction procedure used in 
each case. In all cases, {k) = 60. 



each oscillator's neighbors, we also observe that condi- 
tion (A) leads to the emergence of a bipartite-like net- 
work where low frequency oscillators are mainly coupled 
to high frequency oscillators, as shown in panels (c) and 
(d) for increasing 7. 
We have verified that fulfillment of condition (A) leads 



to a first-order transition for a very wide class of dis- 
tributions of the oscillators' natural frequencies. Fig- 
ure [3]^a)-(b) shows, for instance, the case of the asym- 
metrical Rayleigh distribution, conveniently re-scaled to 

the interval [0,1], given by p{oj) = ■^e~^ (see inset 
Fig. |3ja) with a = 0.25). The results highlight the pres- 
ence of the abrupt transition [Figjsj^a)], together with the 
emergence of a clear frequency-degree correlation, well 
described, again, by the function f{uj) [Figjsj^b)]. 

We now move to discussing several ways to soften con- 
dition (A), while still keeping the abrupt character of the 
transition. For the first extension, we consider again the 
case of a uniform frequency distribution in the interval 
[0,1]. In this case, condition (A) can be relaxed as fol- 
lows: 

(A') for each oscillator z, all nodes j belonging to JV{i) 
satisfy \uji — {ojj)\ > 7c, where (...) indicates the 
average value over the ensemble JV{i). 

The new condition is tantamount to softening condition 
(A) to the local mean field of frequency differences in the 
neighborhood of each network node. Figure [3jc)-(d) re- 
ports the results for networks obtained with a modified 
construction procedure, in which pairs of randomly se- 
lected nodes are now linked if the value of \uji — {ujj)\ 
(averaged over the set of nodes j already linked to node 
i, and the one candidate to be further linked), exceeds a 
gap 7. Again, an explosive transition occurs [Fig. [sj^c)] 
in correspondence to the emergence of frequency-degree 
correlations [Fig. |3jd)]. 

Furthermore, it is worth noticing that a strict appli- 
cation of condition (A) for non uniform frequency distri- 
butions implies that oscillators at different frequencies 
would in general have a different number of available 
neighbors in the network. That's why a natural exten- 
sion of condition (A) is to consider a frequency-dependent 
gap 7(cj) defined by J^^^ p{uj')duj' = Z. Panels (e)-(h) of 
Fig. [3] report the case o? two symmetrical frequency dis- 
tributions, limited to the same frequency range [0, 1], for 
the sake of comparison: i) di Gaussian distribution cen- 

_ (u;-0.5)^ 

tered at a; = 0.5, and given by p{uj) = ^^^ e 2^2 ^ 
with a = 0.13 [Fig. |3je)-(f)], and ii) a bimodal distri- 
bution derived from the same Gaussian, and given by 



p{uj) = 



e 2^2 if^ < 0.5, andp(cj) = 



27T 



Otherwise [with a = 0.23, Fig. [3];g)-(h)]. The gap con- 
dition for the construction of the network is now to fix 
the value of Z, and accept the pairing of nodes when 
\(jUi — ujj\ > ^[j{ijUi) + 7(^j)]- Once again, an explosive 
transition is obtained, with pronounced frequency- degree 
correlation features, as long as p{uj) is symmetrical. 

Finally, we show how our findings allow expressing 
magnetic-like states of synchronization in a such ensem- 
ble of networked oscillators, provided that the coupling 
strength is set inside the hysteresis region of the first- 
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FIG. 4: (Color online) (Left panels) Time evolution of the 
parameter r(t) (see text for definition) for d = 0.004 (top 
panel, outside the hysteresis region), and d = 0.009 (bottom 
panel, inside the hysteresis region). The initial frequencies of 
the oscillators are taken from a uniform distribution in the 
interval [0, 1]. 7 = 0.49 and (k) = 40. Up = 1.0 and Dp = 
0.0005 (bottom red line). Dp = 0.005 (middle blue line), and 
Dp = 0.02 (top black line). The pacemaker is active from t = 
50 to t = 350, as marked by the vertical dashed lines. (Right 
panel) Colormap of aS = (r(t))t>350 (coded as indicated in the 
color bar), showing the region of the parameter space Dp-ujp 
where the magnetic-like state of synchronization is maintained 
after removal of the pacemaker. 

order phase transition. For this purpose, we consider 
again the case of an initial uniform distribution of the 
oscillators' frequencies, and modify system ([T]) as follows: 
^ = Ui^DpS\n{(l)p-(l)i)^dY,f^^ aij sin((/)^- -^^), where 
Dp is the strength of a unidirectional connection to an 
external pacemaker (equal for all oscillators in the net- 
work), and is the phase of the pacemaker obeying 
= ujp. Initially, the system freely evolves in the un- 
synchronized regime, with Dp = 0. The pacemaker is 
then switched on, and Dp is selected so that all oscilla- 
tors are entrained to the pacemaker phase. At a subse- 
quent time, the pacemaker is switched off again, and the 
state of the system is monitored. The left panels of Fig.|4] 
show that setting d outside (inside) the hysteresis region 
produces a final state that relaxes to the original unsyn- 
chronized behavior (that stays permanently in a synchro- 
nized configuration for sufficiently large Dp values). The 
right panel of the same Figure depicts the regions of the 
parameter space Dp — ujp for which these magnetic-like 
states can ultimately be produced. 

In conclusion, we provided a condition for the occur- 
rence of abrupt phase transitions in networks of phase 
oscillators, proved its validity in several cases, and dis- 



cussed several extensions. Our study generalizes previ- 
ous results, and extends the possibility of encountering 
first-order phase transitions to a large variety of net- 
work topologies, as well as to a large variety of frequency 
distribution of the oscillators. This indicates practical 
methods for engineering networks able to display criti- 
cal phenomena, and the emergence of dynamical abrupt 
transitions in their macroscopic states. Furthermore, the 
evidence for the emergence of frequency-degree correla- 
tions in connection with these abrupt transitions, may 
shed light on the mechanisms underlying the relation- 
ship between topology and dynamics in many real- world 
systems. 
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